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ABSTRACT 

Because of the double-helical structure of DNA, in which two strands of complementary 
nucleotides intertwine around each other, a covalently closed DNA molecule with no 
interruptions in either strand can be viewed as two interlocked single-stranded rings. Two 
closed space curves have long been known by mathematicians to exhibit a property called 
the linking number, a topologically invariant integer, expressible as the sum of two other 
quantities, the twist of one of the curves about the other, and the writhing number, or 
writhe, a measure of the chiral distortion from planarity of one of the two closed curves. 
We here derive expressions for the twist of supercoiled DNA and the writhe of a closed 
molecule consistent with the modern view of DNA as a sequence of base-pair steps. 
Structural biologists commonly characterize the spatial disposition of each step in terms 
of six rigid-body parameters, one of which, coincidentally, is also called the twist. Of 
interest is the difference in the mathematical properties between this step-parameter twist 
and the twist of supercoiling associated with a given base-pair step. For example, it turns 
out that the latter twist, unlike the former, is sensitive to certain translational shearing 
distortions of the molecule that are chiral in nature. Thus, by comparing the values for the 
two twists for each step of a high-resolution structure of a protein-DNA complex, the 
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nucleosome considered here, for example, we may be able to determine how the binding 
of various proteins contributes to chiral structural changes of the DNA. 
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INTRODUCTION 

The structure of a DNA molecule is often described as a succession of base pairs, 
each represented as a rectangular plane. 1 A knowledge of the relative locations of origins 
positioned within these planes and the relative orientations of the short and long axes of 
the rectangles allows one to determine for each pair of adjacent base pairs in the molecule 
- a so-called base-pair step - the numerical values of six rigid-body parameters, three 
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translational: shift, slide, and rise, and three angular: tilt, roll, and twist. " 

Some forty years ago mathematicians defined a "twist" which, shortly after its 
introduction, was applied to DNA and used, like the step-parameter twist mentioned 
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above, to characterize the secondary structure of the molecule. " This twist was defined 
as the value of a certain integral involving two continuous space curves. In the 
application to DNA, the structure of which at that time was often depicted in terms of 
space curves, 11 one of the curves was taken to be the axis of the double helix, and the 
other one of the strands winding about this axis. One then went on to compute the twist of 
the DNA, a unitless scalar representing the number of times the strand wound about the 
helical axis. 

Here we are concerned with differences in the properties of these two twists, the 
step-parameter twist, and the twist of the preceding paragraph, which, because of its 
connection with the global shape of the helical axis of a closed DNA molecule, a 
plasmid, for example, we shall refer to as the twist of supercoiling. In the next section we 
begin by reviewing the definition of the twist of supercoiling, and its well-known 
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connection with the writhing number and the linking number. ' Then, after 
characterizing two space curves consistent with today's picture of DNA as a succession 
of discrete rectangular planes, we go on to describe a method for the computation of the 
twist of supercoiling for a single base-pair step. We also point out how easy it is to 
compute the writhe for the case of a closed molecule with an axial curve envisioned as a 
succession of line segments connecting the origins. 
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Of particular interest is the difference in properties between the step-parameter 
twist and the twist of supercoiling. We note that in a relaxed, undeformed configuration 
of a DNA molecule, the two twists are expected to be close in value for all base-pair 
steps. However, we find that although translations of the base pairs leave the step- 
parameter twist unchanged, that is not generally the case for the twist of supercoiling. It, 
instead, is sensitive to translational distortions that are chiral in nature. To illustrate the 
point, we compare the values of the two types of twist for certain base-pair steps in DNA 
wrapped around the core of eight histone proteins in a nucleosome. 14 We find for the 
chosen base-pair steps an instance in which the local protein-induced deformation of the 
DNA structure preserves the value of the step-parameter twist but alters the twist of 
supercoiling. 



THE TWIST OF SUPERCOILING 

As we pointed out above, in much of the early theoretical work describing the 
equilibrium configurations of DNA, the atomistic details of the molecule were ignored, 
and instead, the structure of the molecule was, in fact, described in terms of two space 
curves. The tertiary structure of the molecule was represented, as shown in Fig. 1, by the 
shape of a smooth space curve C mirroring the shape of the axis of the double helix. 




Figure 1 : Schematic representation of DNA. The double helical axis is given by curve C and one of the 
helical strands by curve D. For purposes of calculation of the twist of D about C, D is to be thought of as 
being traced out by the head of a vector £Xl(s c ) everywhere normal to the tangent vector t c (5 c ) . 
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A second curve D, identified with one of the helical strands, was that traced out by the 
head of a vector sd(s c ) , where s is a constant, and d(s c )is a unit vector normal to the 
tangent t c (5 , c )to C at a position having an arc length s c along C. 

The t c (s c ) , d(s c ) pairs associated with two nearby points along C having arc 
lengths s c and s c + ds c are, in general, somewhat rotated with respect to each other. The 
relative rotational orientation of the two pairs before and after a small change in arc 
length ds c is such that there exists a single rotation of the initial pair (at s c ) by a small 
angle about some axis that leads to t c (s c + ds c ) and d(s c + ds c ). A vector d£l is defined 
pointing in the direction of this axis of rotation and having a magnitude equal to the small 
angle of rotation. The change dt c (s c ) and dd(s c ) that the vectors t c (s c ) and d(s c ) 
undergo during this change in arc length ds c is simply the cross product of dSlwith the 
vector itself. That is, if v c (s c ) stands for either t c (s c ) or for d(s c ), 

d\ c (s c ) = dSlx\ c (s c ). (1) 

For the case of the tangent, Eq. (1) leads to the equation 



One of the Frenet-Serret equations, the three equations relating the tangent to the 
principal normal n c (s c ) and the binormal b c (s c )(= t c (s c ) x n c (s c )) , allows us to write 
dt c (s c ) = n c (s c )K c (s c )ds c where ^(^c) i s m e curvature of C. The first term on the right- 
hand-side of Eq. (2) then becomes 



The second term in that equation, the one containing the component of dSl along 
t c (s c ) , is proportional to the twist density. Replacing t c (s c )by d(s c )in Eq. (2) and then 
taking the projection of the resulting expression for d£l along the tangent shows that 



d£l = t c (s c )xdt c (s c ) + (d£l-t c (s c ))t c O c ) . 



(2) 



t c Oc) x dt c (s c ) = b c (s c )tc(s c )ds c . 



(3) 



d£l ■ t c O c ) = (d(s c ) x dd(s c )) • t c (s c ) . 



(4) 



The twist T(D,C) , in units of number of turns, of D about a length / of C is 
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T(D,C) = [j-] C \dn-t c (s c ), (5) 



■'CI 



where s r - s r = I . 

If the curves C, given by r c (s c ) , and D, given by r D (s D ) , are closed, the 
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conventional twist is simply related to two other integrals, ' so-called Gauss integrals, 
the linking number L(CJ)) 

L(D,C)=[ ±)§ W'WW-H a (6) 
Vlx)^ |r D (j D )-r c (s c )| 

and the writhing number W(C), or writhe for short, 

fF(C)= — rff c c , c c V c c \ , c c) 'd s c ds'. (7) 

The linking number is an integer, a topological invariant, equal to the number of 
times that the curve D passes through a surface bounded by C. 15 (In computing this sum 
each pass-through is assigned a value of either +1 or -1 according to a convention 
consistent with the form of the Gauss integral.) This integer remains unchanged for all 
distortions in shape of the curves C and D as long as the curves do not intersect each 
other during the distortions. The writhe, a property of closed curve C alone, is a measure 
of the chiral distortion of the curve from planarity. Fuller pointed out that its value is also 
what one would get by averaging, over all orientations of a plane P, the sum of the signed 
self-crossings occurring in the planar curves resulting from the perpendicular projection 
of Con/ 5 . 10 

The connection between the twist, the writhe, and the linking number mentioned 
above is given by the well-known equation 9 (For another derivation, see Appendix A.): 

L(D,C)=W(C)+T(D,C). (8) 
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THE TWIST OF SUPERCOILING FOR THE MULTISTEP DNA MOLECULE 

As we have mentioned, details of the structure of DNA are now more realistically 
represented, not in terms of smooth space curves but as a sequence of base pairs, each 
pair of adjacent ones enclosing a base-pair step. 1 The base pair is pictured in Fig. 2 as a 
rectangular planar slab containing an origin o from which a triad of mutually orthogonal 
unit vectors emanate, a short axis s , a long axis 1 , and a normal n (= s x 1 ). 




Figure 2: The vectors associated with a base-pair plane: an origin o, and a mutually orthogonal triad of unit 
vectors, the short axis s, the long axis 1, and the normal n. 

We begin the application of the concept of the twist of supercoiling, described 
above, to this structure by imagining that there are simply line segments connecting the 
origins of adjacent base pairs. How then can a smooth curve C be chosen that gives such 
a picture of the DNA, a picture that seems to show a curve with a discontinuous change 
in its tangent at each base pair? That is, for each base pair, the / for example, the 
incoming line segment has a unit tangent we callt (; . 1} . The tangent of the outgoing 
segment ist (;) . Both of these vectors are defined in terms of the origins o ; ._ l5 o ; , and o ;+1 of 
the z' th base pair and the two base pairs adjoining it 



L a-i) 



l (0 



o, - 0, 



o, - 0,. 



° +i - o, 



(9) 



(Note: Subscripts enclosed in parentheses label base-pair steps and those not enclosed in 
parentheses label individual base pairs.) 
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We can envision a limiting process, shown in Fig. 3, in which we start with a 
curve along which the tangent changes smoothly from t (; ._ n to t (i) in the vicinity of base 
pair i as one moves along a circular arc lying in a plane spanned by these two vectors, 
i.e., lying in the plane having as its normal b ; given by 




Figure 3: The passage from a smooth curve with circular segments of curvature K to the entirely linear 
segments connecting the origins of the DNA base pairs. 

Then the arc of the circle is allowed to decrease to zero as its curvature increases to 
infinity, with the tangent vectors t (j ._js and t (;) remaining unchanged. This process is 
repeated at each base pair. Thus, base pair i+l has an incoming tangent t (;) and an 
outgoing tangent t (I+1) determined as in Eqs. (9) from the origins of base pairs i+l and 
i+2. When the tangent becomes t (;) the curve abruptly becomes a line segment of length 
| o. +I - o,. | directed along t (i) . At base pair i+l the curvature abruptly becomes infinitely 
large again but with the tangent still changing smoothly from t (i) to t ()+I) , along a circular 
arc lying in a plane with 

V. = ,1 ( ° X 1 (,+1 V (») 

1 x t (/+i) I 
as its normal. 
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We now can define in greater detail the nature of a single step, the z th . At base pair 
i curve C for the step begins at that point in the circular arc where the tangent, call it t . , 
is midway between t (! ._ 1} and t (0 i.e., 

t,= ,! (, - 1)+ ! ( v 02) 

I t (i-i) + 1(0 I 

and ends at the corresponding point at base pair i+l, i.e., where the tangent is 

~ t (0 + Vi) (13) 
It +t I 

The unit vector dj at the beginning of the step we first take to be in the direction 
of the projection of the long axis 1. on the plane containing o ; perpendicular to t . . The 
unit vector d) +1 at the end of the step points along the projection of on the plane 
containing o ;+1 perpendicular to t m : 
d , = ft-b / )b < +(l < -t < xb < )(t < xb < ) 



' |(l ! .-b ; )b i . + (l ! .-t i .xb,.)(t,.xb,.) 
and (14) 
-b, + i)b i+1 - t i+1 xb j+1 )(t, +1 xb i+I ) 



d* +1 = 



-b, + i)b j+1 +(l <+1 - t <+1 xb w )(t, +1 xb i+I )| 



Let us call a] the angle (in radians) that dj makes with b ; , and «' +1 the angle that 
dj +1 makes with b i+1 , or, more precisely, 
cos a) = b, • dj , sin a) = t . • b ; . x d| 

and (15) 
cos a) +l = b M ■ d| +1 , sin a) +l = t M ■ b M x d) +1 

These vectors and angles are pictured in Fig. 4. 
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Figure 4: The vectors involved in the calculation of the twist of supercoiling of the DNA base -pair step 
bounded by the / th and the (;+l) st planes. Shown in (a), the four base-pair plane origins needed for the 
determination of the three vectors depicted in (b), t ( . _ 1} t (i) and t (i+1) , which, in turn, are needed for 
specifying t ; and t ;+1 . These last two are normal to the two planes seen in (c). Each of these planes 
contains a d vector and a b vector. Two of the angles needed for the twist calculation as given by Eq. (16) 
are denoted in (d), an enlargement of (c). 

To determine the twist of supercoiling of the z th base-pair step we carry out five 
rotations during which dj is converted to dj +1 . All of the vectors involved during this 
process can be thought of as emanating from a single point. For each of the rotations of 
d 1 the rotational axis, d£l, also passing through this point, will either (a) lie along a 
tangent which is not changing its direction, or (b) the axis will coincide with one of the 
two constant b vectors. For a-type rotations, according to Eq. (5), the twist of 
supercoiling is simply the angle of rotation while for 6-type rotations with the axis 
perpendicular to the tangent, the twist is zero. One: We start with the plane containing d| 
and b ; . The normal to this plane is t. . Now rotate dj about t. until the angle between d 1 
and b ; changes from a]to aL where oL. = (a] + / 2. Two: Rotate d 1 about b,until 
the normal to the plane containing these vectors changes from t. to t (;) . At this point the 
plane contains b /+1 as well as b ; . Three: Rotate d 1 about t (;) until d 1 makes an angle of 
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a'^with b. +1 . Four: Rotated 1 about b. +1 until the normal changes from t (i) to t i+i . Five: 
Rotate d 1 about t f+1 until the angle between d 1 and b ;+1 changes from a ( ' ;) to a) +l . The 
d 1 vector has now become dj +1 . There is a nonzero twist associated with the a-type 
rotations one, three, and five, since, for each of these steps, the vector dSl is directed 
along the tangent. Rotations two and four, on the other hand, since the axis of rotation is 
perpendicular to the tangent, are 6-type rotations with zero twist. In the case of rotation 
one and five the twist angle is Aal, 1 2 where Aa!., = a) +l - a]. For rotation three the 
angle of rotation, which we call /3 (i) , has as its cosine and sine: cos f3 (i) = b ; • b /+] and 
sin = t (;) • b ; x b,. +1 . The twist of supercoiling of the base-pair step associated with 
these rotations T\ (in units of number of turns) is thus 
. A«' + /?,.. 

^(0= " • (16) 
in 

The beginning and ending d vectors in the rotations leading to the twist given by 
Eq. (16) are those defined in Eqs. (14) and (15), the unit vectors in the direction of the 
projections of the long axes 1. and onto the planes perpendicular to t,. and t i+1 , 
respectively. 

For an unnicked closed DNA molecule with « B base pairs, Eq. (8) tells us that 

Yjl^LipO-WiC). (17) 

i=i 

If one now were to calculate the T (i) 's for the same five rotations but using as the 
d-vectors those derived from the short axes s, and s i+1 , we would find that for a single 
step the two twists, call them, TL and 73« , would be somewhat different. The sum 

"B 

/T ( s n , however, would have exactly the same value as that obtained before, that given 

1=1 

by the right-hand-side of Eq. (17). We define the twist of supercoiling of the step, T (i) , as 
the average of TL and 73« . Clearly, Eq. (17) is also satisfied for T (i) 's so defined. The 
average of a) and a\ will be denoted as a i . 



12 



In the calculation of the twist of supercoiling for the case of a DNA molecule with 
an open helical axis, for the initial step we take t { to be in the direction of t (1) , and for the 
final step we take t njj in the direction of t (njj _ i;) . 

The step-parameter twist of the base pair step can be derived in a similar way. In 
this case the d vectors are the long and short axes themselves. Carrying out the rotations 
starting with 1 ; and ending with gives the same value for the twist as starting with s ; 
and ending with s ;+1 . There is thus no need to do the averaging indicated in the definition 
of the twist of supercoiling. Figure 5 shows the vectors and angles that play a role in the 
determination of the step-parameter twist. 




Figure 5: The vectors needed for the calculation of the step -parameter twist of the same step shown in the 
previous figure. Here knowledge of the direction of the two normals allows the determination of the single 
vector , which lies in each of the two base-pair planes. Also indicated are the two angles needed for the 
use of Eq. (18) for the twist calculation. 

The analog of t ; and t i+1 are the normals n ; and n /+1 , and b, and b,. +1 become 

the single vector, b'^ = n. x n. +1 / 1 n. x n. +1 1 . The step-parameter twist T% of the step (in 

units of number of turns) thus, is, simply, 

La sp 
2k 



r ( ;; = (18) 



where LaZ is the difference between the angle a^that (or s ;+1 ) makes with and 
the angle af that 1 ; . (or s,. ) makes with b^ . 
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THE WRITHE OF THE CLOSED MULTISTEP DNA MOLECULE 

The shape of the closed axial curve of the DNA molecules we are dealing with, 
consistent with the smooth space curves we have described, is a succession of line 
segments connecting the origins of the base pairs. A method for computing the writhe of 
such a segmented closed curve was elegantly put forth many years ago by Levitt. 16 Here 
we get Levitt's result using a different approach. We first note that for this type of curve, 
the Gauss integral (Eq. (8)) for the writhe for a molecule with n B base pairs, and 
therefore, n B basQ pair steps, takes the form of a sum of contributions w (;)0) of all pairs of 
base-pair steps, 

^(O=I>0-,;) • 09) 

ij 

with 

^-(^ff ^iX^y ^ <2o) 

where 

r w)( s (0' s c/)) = r (o( s (o)- r c/)( s c/)) (21) 

is a vector which points from a point on the / h line segment to a point on the z' th . The 
constant vector X (iJ) is perpendicular to each of the tangents. Its magnitude is thus the 
distance of closest approach of the segments. Terms involving a base-pair step with itself, 
and terms involving pairs of adjacent steps are zero. 

Given Eq. (21), Eq. (20) can be cast into the form 

) (2) f , (;) (2) t m Xt, V X (l 



W dJ) 



n pope) — v^x^) — d d (22) 



2^ ^0)^0) IX^+V,,-!^, 



where 5 (i) (l), ^ (0 (2), * 0) (1), and s {j) (2) denote the arc lengths of the endpoints of the 
segments. 
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The three vectors t (i) , t (j) , and r U J) (s (i) ,s (J) ) determine a dihedral angle 
ju {i j) (s {i) ,s {j) ) we define as follows in terms of its sine and cosine 
r o-j)-( r ( 1 j) xt o-)) x ( t (o xr ov)) 



sin 



r a,j) I ( r (iJ) "t(o x ty)) ^23) 



r (u) x *a) II *(') x r ('j) 



( r 0.7) xt O))-( t (0 Xr 0J)) 

cos / / («-,i) = | — ~[7~ — 

I r (;j) x t y) II t (0 x r v,n I 



It follows from its definition that ju (i j) (s (i) ,s {j) ) is the angle between the vectors normal to 
two planes, the one spanned by r (i J) (s (i) ,s U) ) and t (;) , and the one spanned by 

WVa)) and 

In Appendix B, it is shown that the integrand in the integral appearing in Eq. (22) 
is equal to - S 2 jU (ij) I ds (i) ds (J) so that the contribution to the writhe of base-pair step i and 
base-pair step j, w (iJ) , is simply related to the four dihedral angles associated with the 
endpoints of the steps, i.e., 
' 1 , 

(-// ( , i) (5 (l) (2),5 ( , ) (2)) + M{u . ) (5 (0 (2),5 . ) (1)) 

(24) 



\27T 



We note that w (iJ) has the same sign as t (f) x t {J) ■ X (( . 
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COMPARISON OF THE TWIST OF SUPERCOILING AND THE STEP- 
PARAMETER TWIST 

Like T% of Eq. (18), the twist of supercoiling of a step in the form as expressed 
in Eq. (16) is independent of the direction of propagation along the DNA molecule. But 
whereas does not depend at all on the properties of the two steps adjacent to the i th 
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base-pair steps, T (i) depends on the previous step (7-1), and on the following step (i+l). In 
particular, t,. depends on the origin o ; _j and t ,. + , depends on o,. +1 . For Eq. (17) to be 
satisfied for unnicked closed DNA molecules, the twist of supercoiling must depend on 
the way the tangent vectors are changing. In the case of the step-parameter twist, on the 

other hand, the sum /T^ has no particular significance. 

i=i 

Because in relaxed DNA, the normals of the base pairs and the associated t 
vectors are approximately collinear, for each step we expect the two twists to be close in 
value. When the base pairs are forced to undergo translations of a chiral nature, however, 
unlike the step-parameter twist, which is unaffected by pure translations, we expect the 
twist of supercoiling to change. 

We illustrate this by considering a model structure suggested by an examination 
of a particular three-step section of nucleosomal DNA. Table I shows the twist of 
supercoiling and the step-parameters of three adjacent base-pair steps, 37, 38, and 39, 
derived from the X-ray measurements carried out on this molecule. 14 



Step Shift 


Slide 


Rise 


Tilt 


Roll 


rpsp 




T 














37 -0.18 


-0.18 


3.02 


-0.3 


-0.7 


28.6 


29.6 


38 -0.43 


2.58 


3.25 


-2.3 


-18.4 


50.0 


45.6 


39 0.95 


0.86 


3.47 


2.0 


7.0 


28.5 


28.4 



Table I: Step parameters and the twist of supercoiling of steps 
37-38 of the x-ray crystal structure of a nucleosome core 
particle. The unit of each of the three translational parameters is 
angstroms, and of the angular parameters is degrees. 



Before discussing the implications of Table I, we review the meaning of the five 
step parameters Shift, Slide, Rise, Tilt, and Roll. (An expression for the Twist, here 
called 73? , was derived in a preceding section, Eq. (18).) One, first of all, associates 
with each step a virtual base-pair frame characterized by three mutually orthogonal unit 
vectors, a normal n pointing in the direction of n ; + n. +1 , a long axis I making an angle 
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affs , the average of af and af +x , with , which also lies in this plane, and a short axis 
s given by I x n. Then the shearing translations, shift and slide, are given by the 
projection of o M - o ; on s and on T , respectively, and the rise is given by the projection 
of o ;+1 - o ; on n . Tilt is the angle between the normals of the two base pairs defining the 
step multiplied by the sine of af^ and the roll is the same angle multiplied by the cosine 
of . Each of the six step parameters is associated with a base-pair orientation in Fig. 
6. 



n n 




Slide r ii 




Figure 6: Pictorial definitions of the six rigid-body parameters used to describe the geometry of a base-pair 
step. The short axis, long axis, and normal of the virtual mid-plane of the step are denoted by S , 1 , and 
n . In each of the cases, the image illustrates a positive value of the designated parameter. 
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Note that in Table I step 38 is characterized by a roll of -18°, a slide of 2.6 A, 
and a relatively small tilt and shift. This is consistent with a concerted bending and 

1 7 

shearing deformation underlying the superhelical folding of nucleosomal DNA. The 
twist of supercoiling for each of the three steps is also indicated. For step 38 the twist of 
supercoiling is notably less than the step-parameter twist, unlike the two adjacent steps 
where the two twists are close in value. 

To show how differences in values of the step-parameter twist and the twist of 
supercoiling can arise, the simple model we consider is depicted in Fig. 7. 




31 \ I 



21 * ) 
l i t i 

Figure 7: Construction of a model DNA structure characterized by a chiral deformation. In (a) six equally 
spaced and parallel base-pair planes are shown having their origins lying on a line, (b) depicts the structure 
after the bend described in the text is introduced. The six origins are still coplanar, and the viewing 
direction is chosen to be normal to this plane. A translation of base pairs 4, 5, and 6 as a single unit along 
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the viewing direction, depending on the direction of the motion, results either in a structure with a right- 
handed jog ((c)) or one with a left-handed jog ((d)). 

At the start, in Fig. 7(a), there are six base pairs with their origins collinear and 
equally spaced. In addition, the base-pair planes are oriented so that their normals 
coincide with the t-vectors. For each of the five steps the rise is simply equal to the 
distance between the origins, and the step-parameter twist and the twist of supercoiling 
are equal in value. Also, for each of the steps, aside from the rise and the step-parameter 
twist, the other step parameters are zero. The line of base pairs 4,5, and 6 is now rotated 
by some angle about an axis lying in base pair 3 and passing through the origin of that 
base pair, and then by the same angle about an axis lying in base pair 4 and passing 
through the origin of that base pair. The effect of the rotation is to change the values of 
the roll and tilt of step 3. All of the other step parameters are left unchanged. Neither do 
any of the twists of supercoiling of the steps change. After this bend is introduced, the 
origins of the base pairs are no longer collinear, but they all do remain coplanar. Figure 
7(b) shows the new structure with the viewing direction chosen to be normal to the plane 
on which the origins lie. Finally, a translation of base pairs 4, 5, and 6 as a unit in each of 
the two directions perpendicular to the plane on which their origins had been lying gives 
rise to the two chiral structures, a structure with a right-handed jog shown in Fig. 7(c), 
and one with a left-handed jog shown in Fig. 7(d). The final shearing motion alters just 
the slide and the shift of step 3 in the case of the step parameters, and, in fact, the ratio of 
the roll to the tilt equals the ratio of the slide to the shift. But also accompanying this 
introduction of structural chirality is a change in the twists of supercoiling of steps 2, 3, 
and 4. It is this that we now examine in some more detail. 

According to Eq. (16), we can write for the twist of supercoiling for each of the 
three central steps of Fig. 7(c) or 7(d) 
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«5-«2+ZA< 



with i = 2, 3, 4. The sum of these twists of the steps in question is, therefore, 

4 

2X = (26) 

Since for the final structures of this model all of the vectors are unidirectional and the 
vectors b 2 and b 5 coincide with bf and bf , it follows that 

4 r/ sp - rv sp 

YT- = ^-^- (27) 

«=2 ^ 

and that a 2 = af and « 5 = af . Thus 



2X,-Et;;=<=!— . (28) 

i=2 i=2 Z;T 

That is, for this model, the difference in the sums of twist of supercoiling and the 
step-parameter twist for the three central steps depends only on the rise, roll, and slide. In 
particular, the difference does not depend on the values of the step -parameter twists 
originally assigned in the linear structure of Fig. 7(a). It is also found that the values of 
the differences T (i) - for each of the three steps are only mildly dependent on the step- 
parameter twists of the steps. Furthermore, the value of the twist difference when i = 2 or 
4 is much smaller than it is for the case of the middle step, i = 3. With a roll, slide, and 
rise of step 3 of -18°, 2.6 A, and 3.4 A, respectively, and with zero tilt and shift, the 
model gives, for the case of the case of the left-handed jog, in units of degrees 

4 4 

IX) - Z r (o = -6-10°. Also, when T* = 29°, T* = 50°, and 7$ = 29°, we find that 
T {2) ~ T ( i = -0.01°, T (3) - = -6.08°, and T {4) - T* = -0.01°. 
The model structure with the left-handed jog, Fig. 7(d), is seen to have properties 
qualitatively similar to those shown by the nucleosomal steps 37, 38, and 39. The 

39 39 

corresponding twist differences for the nucleosome are ^ T (i) - ^ T$ = -3.1° , 

!=37 i=37 

T — T sp =0 9° T - T sp =-4 4° and T - T sp =-0 2° 

J (37) J (37) ' J (38) J (38) ^t-t ) J (39) J (39) • 
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Appendix A 

We begin the derivation of Eq. (8) by expressing the connection between curve D 
and curve C, by writing for r D , 

r»(s D ) = r c (s' c ) + £d(s' c ), (A-l) 

where s is a constant and d(s' c ) is a unit vector perpendicular to the tangent t c (s' c ) to 
curve C, i.e., d(s' c ) ■ t c (s' c ) = 0. Equation (A-l) implies that the unit tangent t D (s D ) to 
curve D is of the form 



t r (s') + £ 



(A-2) 



where k is the reciprocal of the magnitude of the vector contained in the parentheses. 
These two expressions are substituted into the integrand appearing in the Gauss integral 
form for the linking number, Eq. (6). We then allow s to approach zero. Beyond a certain 
point in this limiting process curve C and curve D no longer can intersect each other, and 
thereafter the linking number remains unchanged. 10 The factor k approaches one as s 
approaches zero. We also find that, in this limit, the terms in the integrand that are linear 
in s are zero. Thus one can write for the linking number 



L(D,C) = W(C) + lim 



d(, c )x^)-t c (^) 

§- C , 2 ^n ds' c ds c (A-3) 

W (|r c (5 c )-r c (4)| 2 +s 2 ) 



£ 



where W(C) is the writhing number as given in Eq. (7). Furthermore, because 

s s 
Hm = lim 

~ (jr c (, c )-r c (^)| 2 + ^f ^°((, c -4) 2 + ^ 2 J (A-4) 

= 2S(s c -s' c ), 

where S(s c - s' c ) is a Dirac delta function, the integration over s' c in the second term of 
Eq. (A-3) can be readily carried out: 
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lim 

£->0 



Iff 



d(s c )x dd ^l.t c (s' c ) 



d,S r 



(|r c (, c )-r c (4)| 2 +^ 2 ) 



3/2 ""C^C 



f 1 A 



V 2^y 



$ 8{s c - 4) d(* c ) x ^1 • t c (4) (A-5) 



\2x ) 



|d(5 c )x^l.t c (j c )d! c . 



This completes the proof of Eq. (8). 



Appendix B 

We begin by deriving an expression for - dju (i j) I ds ( j) . Since, for example, 
Scot// _ d cot ju dju 



8s, 



U) 



d ju 8s u 



(B-l) 



we see that 
dju 



8s (J) 



. 2 0COt« 

sin ju — 

8s U) 



(B-2) 



Given the definition of the dihedral angle (Eqs. (23)), and the fact that 
t (i) x t 0) • r (iJ) is independent of s {i) and s {J) , this last equation can be rewritten as 

8ju I r ( ,-j) I 2 (t (0 xt (j) ••"(«,,)) 



ds (j) 



r o-j) r i, t (,-) xt o-) • r Q-, 7 ), 

\ r (iJ) Xt U) I l t (0 Xr (i,y) I 



5 



as. 



(./) 



(B-3) 



After noting that the explicit dependence of r (;j) on s (i) and s {j) as given in 
Eq. (21) leads to the fact that 
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r (',7) I ~l X (< J) I +S U) +S U) 2 ( t 0) '^^(O^O) 



t (o xt y) ' r oj) _t (o xt o) ' x oj) 



^x^l'HX^r+lt^xt^l 2 ^, (B-4) 



I r a,j) xt u) I 2- ' X ('j) I 2 + 1 *(«•) xt u) I 2 S l) 

( r (U) xt oo)-( t (o xr (u)) = -(( t (o- t o-))l X (u) I 2 +l t (o xt (v) I 2 Wo-)) 

We find, after performing the indicated differentiation with respect to s (j) in 
Eq. (B-3) and simplifying the resulting expression, that 

d/j = (t( f ) xt -) -Xq-j)) (^(o -(*(,-) • t o)>o)) 

& w (l^l'+l^xto-)! 2 ^))!^-)! • (B " 5) 



Differentiating Eq. (B-5) again, this time with respect to s (i) , yields 

9 V(,,y) _ t (o xt (./)' X o-j) 



(B-6) 
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